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ABSTRACT 

We report on a new survey for z ~ 4.5 Lyman-a sources, the Large Area Lyman 
Alpha (LALA) survey. Our survey achieves an unprecedented combination of volume 
and sensitivity by using narrow-band filters on the new 8192 2 pixel CCD Mosaic 
Camera at the 4 meter Mayall telescope of Kitt Peak National Observatory. 

Well-detected sources with flux and equivalent width matching known 
high redshift Lyman-a galaxies (i.e., observed equivalent width EW > 80A, 
2.6 < (line + continuum nux)/(lCr 17 erg cm" 2 s" 1 ) < 5.2, and <5(EW)/EW < 0.25) 
have an observed surface density corresponding to 11000 ± 700 per square degree 
per unit redshift at z = 4.5. Spatial variation in this surface density is apparent on 
comparison between counts in 6561 ± 40A and 6730 ± 40A filters. 

Early spectroscopic followup results from the Keck telescope included three sources 
meeting our criteria for good Lyman-a candidates. Of these, one is confirmed as a 
z = 4.52 source, while another remains consistent with either z = 4.55 or z = 0.81. We 
infer that 30 to 50% of our good candidates are bona fide Lyman-a emitters, implying 
a net density of ~ 4000 Lyman-a emitters per square degree per unit redshift. 
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1. Introduction 

More than three decades ago Partridge and Peebles (1967) predicted that galaxies in their first 
throes of star-formation should be strong emitters in the Lyman-a line. Their predictions were 
optimistic, based on converting roughly 2% of gas into stars in 3 x 10 7 years in Milky Way sized 
galaxies, which translates into a luminosity of 6 x 10 44 ergs _1 . These objects are also expected to 
be common - if all the L* galaxies have undergone a phase of rapid star-formation one should see a 
surface density of about 3 x 10 3 x (At/ (3 x 10 7 yr))deg~ 2 (Pritchet 1994). Searches based on these 
expectations did not detect Lyman-a emitters (LAEs). (See review by Pritchet 1994; Koo & Kron 
1980; Pritchet & Hartwick 1987, 1990; Cowie 1988; Rhea et al 1989; Smith et al 1989; Lowenthal 
et al 1990; Wolfe et al 1992; De Propris et al 1993; Macchetto et al 1993; M0ller & Warren 1993; 
Djorgovski & Thompson 1992; Djorgovski, Thompson, & Smith 1993; Thompson, Djorgovski, & 
Trauger 1992; Thompson et al 1993; Thompson, Djorgovski, & Beckwith 1994; Thommes et al 
1998.) 

Only recently have Lyman-a emitters been observed, albeit at luminosity levels roughly a 
hundred times lower than the original prediction. These Lyman-a emitters have been found from 
both deep narrow band imaging surveys (Cowie & Hu 1998; Hu, Cowie & McMahon 1998; Hu, 
McMahon, & Cowie 1999; Kudritzki et al 2000), and from deep spectroscopic surveys (Dey et 
al 1998; Manning et al 2000; but see Stern et al 2000). Weak Lyman-a emitters have also been 
found through targeted spectroscopy of Lyman break objects (e.g., Steidel et al 1996, Lowenthal 
et al 1997). The lower luminosity in the Lyman-a line may be because of attenuation by dust if 
chemical enrichment is prompt; or because the star-forming phase is more protracted; or because 
the star-formation happens in smaller units which later merge. The first two scenarios will give a 
smaller equivalent width than early predictions, while the last scenario results in low luminosities 
but high equivalent width. 

Dust effects are expected to be severe — even a small amount of dust can greatly attenuate 
this line, because it is resonantly scattered. However, two factors can help the Lyman-a photons 
escape. If Lyman-a photons are produced in diffuse regions of a clumpy interstellar medium, they 
can simply scatter off the dense clumps and escape (Neufeld 1991), and some geometries can even 
lead to an increase in the equivalent width of the line. Secondly, energetic winds are seen in low-z 
Lyman-a emitters (Kunth et al 1998). These can displace the neutral gas and doppler-shift the 
peak wavelength of the resonant scattering, thereby reducing the amount of scattering and the 
path length for interaction with dust. 

Detailed predictions for luminosities and surface densities of LAEs using a Press-Schechter 
formalism and exploring a range of dust obscuration and star-formation time scales have been 
explored by Haiman and Spaans (1999), who are able to reproduce the surface densities of LAEs 
reported by Hu et al (1998) with a wide range of models. In order to narrow down the range of 
possibilities and characterize the high redshift Lyman-a population, better statistics over a wide 
range of flux and source density are needed. 
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2. Narrow Band Imaging Survey 

An efficient search for Lyman-a emitters (and other emission line galaxies) was started in 
1998 using the CCD Mosaic Camera at the Kitt Peak National Observatory's 4m Mayall telescope. 
The Mosaic camera has eight 2048 x 4096 chips in a 4 x 2 array comprising a 36' x 36' field of view. 
The final area covered by the LALA survey is 0.72 square-degrees in two MOSAIC fields centered 
at 14:25:57 +35:32 (2000.0) and 02:05:20 -04:55 (2000.0). Five overlapping narrow band filters of 
width FWHM« 80Aare used. The central wavelengths are 6559, 6611, 6650, 6692, and 6730A, 
giving a total redshift coverage 4.37 < z < 4.57. This translates into surveyed comoving volume 
of 8.5 x 10 5 comoving Mpc 3 per field for H = 70kms- 1 Mpc -1 , tt = 0.2, A = 0. About 70% of 
the imaging at z ~ 4.5 is complete, and an extension of the survey to z > 5 is planned. In about 
6 hours per filter per field we are able to achieve line detections of about 2 x 10~ 17 ergcm~ 2 s _1 . 
The survey sensitivity varies with seeing. Broadband images of these fields in a custom B w filter 
(A = 4135A, FWHM = 1278A) and the Johnson-Cousins R, I, and K bands are being taken as 
part of the NOAO Deep Widefield Survey (Jannuzi & Dey 1999). 

The images were reduced using the MSCRED package (Valdes &; Tody 1998; Valdes 1998) 
in the IRAF environment (Tody 1986, 1993), together with assorted custom IRAF scripts. In 
addition to standard CCD data reduction steps (overscan subtraction, bias frame subtraction, 
and flatfielding) , it was necessary to remove crosstalk between pairs of CCDs sharing readout 
electronics and to remove a ghost image of the telescope pupil. Astrometry of USNO-A catalog 
stars was used to interpolate all chips and exposures onto a single coordinate grid prior to 
combining final images. Cosmic ray rejection is of particular importance in a narrowband search 
for emission line objects. We therefore identified cosmic rays in individual images using a digital 
filtering method (Rhoads 2000) and additionally applied a sigma clipping algorithm at the image 
stacking stage. Weights for each exposure were determined using the method of Fischer & 
Kochanski (1994), which accounts for sky level, transparency, and seeing to optimize the signal to 
noise level for compact sources in the final image. 

Catalogs were generated using the SExtractor package (Bertin & Arnouts 1996). Fluxes were 
measured in 2.32" (9 pixel) diameter apertures, and colors were obtained using matched 2.32" 
apertures in registered images that had been convolved to yield matched point spread functions. 

3. Spectroscopic Observations 

Spectroscopic followup of a cross-section of emission line candidates was obtained with the 
LRIS instrument (Oke et al 1995) at the Keck 10m telescope on 1999 June 13 (UT). Two dithered 
exposures of 1800 seconds each were obtained through a single multislit mask in good weather. 

These data were reduced using a combination of standard IRAF ONEDSPEC and 
TWODSPEC tasks together with the homegrown slitmask reduction IRAF task "BOGUS" (Stern, 
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Bunker, & Stanford, private communication) for reducing LRIS data. 

4. The emission line source population 

Our imaging data yield numbers of sources as a function of their fluxes in several filters, 
from which we can construct number densities as a function of magnitudes, colors, and equivalent 
widths. In order to gracefully handle sources that are not detected in all filters, we have chosen to 
use "asinh magnitudes" (Lupton, Gunn, & Szalay 1999), which are a logarithmic function of flux 
for well detected sources but approach a linear function of flux for weak or undetected sources. 
Figure |l] shows the color-magnitude diagram for the 6559A (±40A) and R filters. The color scatter 
achieved for bright sources (R < 22) is 0.10 magnitudes (semi-interquartile range). This includes 
the true scatter in object colors, and is therefore a firm upper limit on the scatter introduced by 
any residual systematic error sources, which we expect to be a few percent at worst. 

To sharpen our focus on the high redshift Lyman-a population, we identify the range of 
parameter space occupied by known z > 3 Lyman-a emitters. These sources have typical observed 
equivalent widths ;> 80A and line+continuum fluxes ^5x 10~ 17 erg cm -2 s _1 (adjusted to z = 4.5 
and measured in an 80 A filter). 

We further want to restrict attention to sources with sufficiently reliable detections that the 
number of false emission line candidates is a small fraction of the total candidate sample. The 
data set includes many continuum sources (~ 1.4 x 10 4 in the flux range above). It is expected to 
contain a few hundred emission line sources, based on earlier source counts from smaller samples 
(Hu et al 1998) and on our own findings (below). The fraction of continuum sources with measured 
equivalent widths above some threshold EWo can be calculated as a function of signal to noise 
level n and filter width AA: The measured equivalent width for a source with no line emission 
or absorption in the narrowband filter will be ± AA/n. Thus, a false positive becomes an ma 
event, where m = n x EWo/AA. Our sample would be reasonably safe from contamination with 
m = 3, which would correspond to about 20 false positives in a sample of 14000 sources; and 
very safe for m ^ 4, corresponding to < 1 false positive in 14000 sources. We have conservatively 
chosen detection thresholds n = 5 and EWo = AA = 80A, which gives m = 5. This keeps the 
number of false positives small even when the foregoing analysis is expanded to include errors 
in the continuum flux (which increase photometric error in the color by a factor ~ \/2) and a 
realistic distribution of equivalent widths for the \ow-z galaxy populations. As a further check on 
our sample, we use the photometric color error estimates from SExtractor to demand that the 
source be an emission line source at the 4cr level. 

Combining all of these requirements, the final criteria for good candidates in our survey 
become EW > 80A, 5(EW)/EW < 0.25, and 2.6 < f 17 < 5.2, where f 17 = f/(W~ 17 erg cm" 2 s" 1 ). 
There are 225 such sources detected in the 6559A filter in a solid angle of 0.31 square degree and 
redshift range 5z = 0.07. This corresponds to 11000 good candidate Lyman-a emitters per square 
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degree per unit redshift. The precise upper flux cutoff used is somewhat arbitrary, but including 
sources with larger fluxes in our list of good candidates has little effect on the total source counts. 

The final 6730A image had a broader PSF than the 6559A image, and consequently a 
reduced effective sensitivity. In addition, it samples a slightly higher redshift, resulting in an 
0.08 magnitude increase in distance modulus (qo = 0.1, A = 0). Accounting for both effects, the 
matched flux range is 3.45 < fn < 5.2 at 6559A and 3.2 < f 17 < 4.8 at 6730A. The corresponding 
counts are 70 candidates at 6730A, and 104 candidates at 6559A. Thus, the source density in 
directly comparable luminosity bins varies by a factor of 1.5, a difference that is significant at 
about the 3a level. We interpret this as a likely signature of large scale structure in the Lyman-a 
emitter distribution at z ~ 4.5. The comoving distance between the centers of the two redshift 
slices is 72Mpc, while the thickness of each slice is 36Mpc and the transverse comoving size 89Mpc 
(again assuming Hq = 70 kms -1 Mpc -1 , Q = 0.2, A = 0). Comparable variations have been 
observed in the density of Lyman break galaxies at z « 3.1 (Steidel et al 1998). 

5. Spectroscopic results 

Our spectroscopic followup sampled a wide range of flux and equivalent width in order to 
characterize the different populations of sources with extreme narrowband colors, and to thereby 
tune candidate selection criteria for future spectroscopic observations. Those sources with 5a 
narrowband detections and photometrically measured EW > 65A were all confirmed as emission 
line objects at the narrowband wavelength. In total, we detected two [O III] A5007 emitters at 
z = 0.34, three [O II] A3727 emitters at z = 0.77 and z = 0.81, one confirmed Lyman-a (1215A) 
emitter at z = 4.516, and one source that could either be Lyman-a at z = 4.55 or [O II] A3727at 
z = 0.81. We also found a z = 2.57 galaxy with emission lines of C IV A 1549, He II A 1640, 
and O III A1663. This source may have been included in the narrowband candidate list (with 
EW = 36A) because of a weak spectral break around 6560A. Finally, we found a few serendipitous 
emission line sources in the slit spectra. One of these is a single-line source with large equivalent 
width, possibly Lyman-a at z = 3.99. Spectra of two interesting sources are shown in figure ||. 

6. Discussion: The Lyman-a source population 

By combining our imaging survey with these spectroscopic results, we can estimate the source 
density of Lyman-a emitters passing our selection cut. Our spectra included three sources fulfilling 
the criteria given above for good candidates. Of these, one was confirmed as a z = 4.52 Lyman-a 
source. A second remains a candidate z = 4.55 source, but is more conservatively interpreted as a 
z = 0.81 [O II] emitter on the basis of a rather strong continuum on the blue side of the line. The 
third is a clear z = 0.34 [O III] emitter. We therefore estimate that roughly 1/3 to 1/2 of the good 
candidates will be confirmed as Lyman-a sources, yielding ~ 4000 emitters per square degree per 
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unit redshift. This is compatible with earlier measurements from smaller volumes (Hu et al 1998) 
after accounting for differences in flux threshold. 

Our measurement is distinct from previous efforts in the field for its basis in a large number of 
candidate emitters. Poisson errors in our source counts are of order ±7%. This is smaller than the 
variations observed in the comparison of two filters (of order ±40%). By combining observations 
in multiple fields, we will be able to average over local fluctuations in number densities effectively. 
When completed, the LALA survey will yield comoving volume of ~ 2 x 10 6 Mpc 3 (§|) and 
a sample of several hundred LAEs, and will allow the luminosity function, equivalent width 
distribution, and correlation function of this population to be determined for the first time. 

We thank Andy Bunker and Steve Dawson for help with the spectroscopic observations and 
Frank Valdes for writing and helping with the MSCRED package in IRAF. JER's research is 
supported by a Kitt Peak Postdoctoral Fellowship and by an STScI Institute Fellowship. SM's 
research is supported by NASA through Hubble Fellowship grant # HF-01111.01-98A from 
the Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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Fig. 1. — A sample LALA survey color- magnitude diagram, for narrowband 6559A mag nitude (80A 
bandpass) and 6559A — R color. The zero point of the narrowband filter is calibrated using R band 
standard star fluxes (Landolt 1992), so that m(6559A) = R for a pure continuum source of average 
spectral slope. The box in the lower right hand corner shows the color and magnitude constraints 
for good Lyman-a candidates, as described in the text. Additional constraints on the color error 
eliminate some objects meeting the color-magnitude cut; the final sample of good candidates for 
this filter is marked by the heavy points. 
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Fig. 2. — Keck spectra of two LALA emission line candidates. Top: The largest equivalent width 
source in the spectroscopic sample, a z = 0.34 OIII emitter with observed EW ^ 1400A. This 
source has eleven detected emission lines, and an [O III] A5007A line luminosity of 2.5 x 10 41 ergs -1 . 
It demonstrates our ability to find unusual objects through our large volume coverage. Bottom: A 
confirmed z = 4.516 Lyman-a source. This object has a line flux of 1.7 x 10~ 17 erg cm -2 s _1 and 
an equivalent width of 84A. The line is asymmetric and has a strong continuum decrement from 
the red to the blue side, both of which are expected for high-redshift Lyman-a emitters (Stern & 
Spinrad 1999). It demonstrates our ability to find the faint Lyman-a emitter population. 



